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1. Introduction 

Bhabha scattering has been measured with great accuracy by the experiments at 
the LEP and SLC coUiders. This has led to many works by numerous authors being 
devoted to calculating the theoretical cross section for this process as accurately as 
possible [1, 2]. Yet, the current accuracy of theoretical predictions for large angle 
Bhabha scattering at LEP2 is still not completely satisfactory [3], and limits the 
bounds that can be set on some kinds of potential new physics effects, for instance, 
those coming from large extra dimensions [4]. 

Future linear electron positron colliders, e.g. the proposed TESLA accelerator, 
will not have monochromatic beams because the electrons and positrons can emit 
beamstrahlung before they collide. Measurement of the resulting acoUinearity angle 
in large angle Bhabha events provides a way to determine the luminosity spectrum [5, 
6]. 

Large angle Bhabha scattering is also important at electron-positron collid- 
ers running at centre of mass energies of a few GeV, such as BEPC, VEPP-2M, 
DAPHNE, and the B-factories PEP-11 and KEK-B, where it is used to measure the 
integrated luminosity [7]. At present, the theoretical uncertainty on the differential 
cross section of this process is one of the hmiting factors on the precision of the 
luminosity determination [8]. 

In order to improve the theoretical predictions, it is necessary to include higher 
order radiative corrections. In this paper, we consider only pure QED corrections. 
These fall in different orders of magnitude aL, a, a^LP', a'^L, a^L^, where a is 
the fine-structure constant and L = log{s/ml). The large logarithm L is related 
to collinear divergences that would appear if the electron mass rUf, were zero. The 
O (a^L) terms are so far only partially known. Arbuzov, Kuraev and Shaikhatdenov 
have calculated the contributions from soft one- and two-photon bremsstrahlung, 
squared one-loop graphs, and the interference between two-loop vertex graphs and 
tree level terms [9]. However, they did not calculate the contribution due to two- loop 
box graphs, which has to be included to complete the O (a'^L) correction. 

In principle, it should be possible to extract the missing two-loop box contribu- 
tion from the work of Bern, Dixon and Ghinculov [10] , who have presented a complete 
formula for the two-loop virtual corrections to Bhabha scattering. However, these 
authors set the electron mass to zero and regularized all infrared and collinear di- 
vergences dimcnsionally, with d = 4 — 2e, unlike the authors of ref [9], who used the 
traditional method of an electron and a photon mass. For this reason, the results 
obtained by the two groups cannot be directly combined. 

If one is only interested in the logarithmically enhanced term, one does not 
actually need the full result of ref. [10], but only the terms containing powers of 1/e. 
The correspondence between poles in e of the dimensional regularization scheme and 
logarithms log(m-y) and log(me) of the massive regulator scheme has been worked 



1 



out in detail at the one-loop level [11]. However, it is not a priori clear how to extend 
this to two loops. Therefore, in this paper, we take a slightly different approach. 
We rely on an elegant classification, proposed by Catani, of the singularities of on- 
shell two-loop amphtudes in QCD [12]. His results were confirmed for the QED 
case by [10]. Catani's formalism was developed using dimensional regularization. 
However, it ought to be valid in both regularization schemes. We shall translate this 
formalism into the massive regulator scheme, and then use it to predict the O (a^L) 
contribution to the cross section for large angle Bhabha scattering. We limit ourselves 
to graphs without photon vacuum polarization insertions, since those can be treated 
as a separate class [13]. 

Our result should be directly applicable to the energy range of a few GeV, 
where non-QED effects, from graphs involving Z-boson exchange, are still very small. 
There are several Monte Carlo programs that were specifically designed for this en- 
ergy range. These include the programs BABAYAGA and LABSPV [7] and LAB- 
SMC [14]. (Monte Carlo programs written for LEPl and LEP2 are reviewed in 
refs. [2, 3]). The precision of these programs depends on the energy and on the 
details of the cuts applied. In the case of LABSMC, under typical conditions, the 
uncertainty is estimated by the authors to be around 0.2% [14], mainly due to missing 
terms of order O (a^L). The result obtained in this paper will allow this uncertainty 
to be reduced. 

2. Infrared factorisation formulae 

A generic (renormalized) QED matrix element can be expanded as a series in a as 
follows, 

where represents the i-loop contribution and where the overall power n may 

be half-integer. According to Catani [12], these amplitudes obey the following fac- 
torization formulae, 

= /(2)|>j(o)) + + (2.2) 

where and are operators that depend on the regularization scheme and 
contain all of the singularities of the infrared regulator. The remainders \Ai^^^'-^^"') 
and \A4^'^^'f"^) are finite. In dimensional regularization, for QED processes with 
massless electrons. 
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where \ij = —1 Hi and j are both incoming or outgoing partons and \ij = 
otherwise, and Cj is the electric charge (minus the electric charge) of an outgoing 
(incoming) radiating particle with momentum pi. Eq. (2.3) is obtained from the 
QCD result [12] with the replacements Ca — > 0, Ci;' — > 1, — > 1 and T^-Tj — > ejCj. 

Because we exclude graphs with photon vacuum polarization insertions, we ne- 
glect terms proportional to the (5 function. The operator J^^^ is then given by 

^ _lj(i)j(i)+jj(2)^ (2.4) 
where H^'^^ contains only single poles and has the form 

where 

= ^ - 3C2 + 6C3. (2.6) 

Ref. [12] does not give a general formula for H^'^\ and at present it is necessary to 
derive it from an explicit two-loop calculation, such as the two-loop QCD computa- 
tion of the electromagnetic form factor of the quark [15]. Wc note that in Ref. [12] 
the summation over radiating pairs i and j is not explicitly included in the definition 
of H^'^\ However, it has been found to reproduce correctly the results of explicit 
calculations in QED [10] and QCD [16]. We also note that the terms of C (1) in if^^^ 
are presently a matter of choice and can be altered by a redefinition of | Al^^^'-'''"). 



3. The massive photon and massive electron regularisation 
scheme 

The factorization formulae (2.2) should hold for any choice of infrared regulator. 
Therefore, we can find a translation between schemes by comparing explicit calcula- 
tions. For example, the one-loop correction to the electron photon vertex (see, e.g. 
Ref. [17]), serves to fix I^^^ in the scheme where the photon has mass A and the 
electron has mass m. We write the Dirac form factor Fi[q^) of the electron photon 
vertex 

= l,F,{q^) + ^^M-^''^2(g') (3.1) 

as-*- 

F,{q^)^l + ^F['\q^) + (^^)'F['\q^)+0(a') . (3.2) 

^In order to follow more closely the notation of Ref. [9], we expand in a/vr rather than Q;/(27r) 
in the remainder of this paper. Therefore, the definitions of J^^^ and H^'^^ differ by factors of 1/2 
and 1/4 respectively from those in the previous section. 
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By equating 

Fi'\q') = I^'\q')+F}'^'^'-{q'), (3.3) 

we see that 

I''' (g^) = (L - 1) + l) - - + g , (3.4) 

where Lx — log(A^/m^), L — log(— g^/m^) if g'^ < and L — log(g^/m^) — in if 
> 0. There is a possible ambiguity in assigning the constant pieces to I^^^ or 
p{i),fzn^ Our choice corresponds to F^;^^'^^"' — 0. Similarly, by examining the two- 
loop vertex, we can fix H^'^\ To do this, wc take the large scale limit of the O {e'^) 
two-loop vertex correction F^\q^) computed by Barbieri, Mignaco and Remiddi [17] 
(without the contribution from the vacuum polarization graph). Up to terms that 
do not depend on the regulators, we expect that, 

F?\q') = \l^'\q')I^'\q') + H^'\q') + O {!) . (3.5) 

We find that H^'^\q^) is proportional to a single large logarithm 

Jf(^)(g2) = iLi/(2) (3.6) 

where H^"^^ is given by Eq. (2.6). In fact, one might wonder whether changing the 
constants in I^^\q^) would give rise to a different H^'^\ This is not the case since 
any change is absorbed by the necessary alteration to F^^^'^'^'", 

Armed with these operators using mass regularization, we can compute the large 
logarithmic corrections to the two-loop contribution to Bhabha scattering keeping 
the electron mass and using a small photon mass as the infrared regulator. For 2 — > 2 
scattering, there are six radiating pairs, so that 

/(I) = 2 (/(^)(s) + I^'\t) - I^'\u)) , (3.7) 

and 

H(2) ^ 2 (Jf(2)(s) + Jf(2)(i) - H^^){u)) . (3.8) 

The second order virtual contribution to Bhabha scattering comes from the 
square of the one-loop graphs and the interference of tree and two-loop graphs. We 
can write this contribution in factorized form as 

' ] Ayy, (3.9) 



d(7o Vtt/ 

where the lowest order differential cross section is given by 



dao = y (^-^^^Xdn, (3.10) 
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with X = —t/s and 1 — x = —u/s. Up to constant terms, 

with J^^^ and i?^^^ given by Eqs. (3.7) and (3.8) respectively. The second term 
involves the one-loop virtual contribution 5v (see, e.g. [18, 19] and references therein). 
In the notation of ref . [9] , it is given by. 



m 



4 log y ( 1 - L + log 



1 — x 



X 



L2 + 2Llog 



1 — X 



X 



log^(a;) 



+ log'(l-x) + 3L-4 + /(x), 



(3.12) 



with^ 



f{x) = {l-x + x') 



2\-2 



TT 

12' 



- 8x + 27x^ - 26x^ + 16x^) 



+^(-2 + 5x - 7x2 ^ 5^3 _ 22;4-) iog2(2 - x) + ^a;(3 -x-2>x^^ Ax^) \og^{x) 

+ -(6 - 8a; + 9x^ - 3a;^) log(a;) - -a;(l + x^) log(l - x) 
2 2 



+-(4 - 8x + 7^2 - 2x^) log(a;) log(l - x) 



(3.13) 



Expanding Ayy yields all terms containing at least one power of the large logarithm 
L = log(s/m2) as well as all logarithms of the photon mass regulator. We have 
checked that in the small angle limit, a; ^ 0, to the logarithmic accuracy we are 
working at, Eq. (3.11) reduces to the known result [20] 



Avv^6{Fi'\t)y + AFi'\t), 



(3.14) 



which follows from a generalized eikonal representation [21] for the Bhabha scattering 
amplitude for small angles. 



4. The scattering cross-section at O (a'^^ 

To make a physical prediction, the double virtual contribution must be combined 
with the one-loop contribution with single soft emission and the tree level double 
soft emission. The second order correction to the one-loop virtual photon emission 
corrected cross section, due to the emission of a single real soft photon having energy 
less than As, can be written down in the factorized form [9] 

^-Ss-Sv=i-) Asv, (4.1) 



d(To TT TT VtT/ 



^There are some misprints in the formula for f{x) in ref. [9]. 
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where 5y is given in Eq. (3.12) and 

'niAe\ ( ^ ^ , fx 



- 41og (l - 1 + log (^)) + + 2Llog (j^) + log^(.) 

27r2 

- log2(l - x) - — + 2 Li2(l -x)-2 Li2(a;). (4.2) 



Here, £ = \/'5/4 is the energy of the electron and positron beams. The dilogarithm 
function is defined as 

dy 

y 



Li2(x) = -/^log(l-y). (4.3) 



The contribution from two independently emitted soft photons each with energy 
1^1, 1^2 < is given by [9] 

da^^ 1 /aX 



{^)'sl.{^)\ss. (4.4) 



d(To 2 UJ ^' 

where the statistical factor 1/2! is due to the identity of photons. 

The combination Ass + ^sv is written in expanded form in Rcf. [9] where terms 
like L*^, L^Lx, L^, L^Lx, L'^L\, LLx and LL\ arc produced. All of these terms 
precisely cancel against similar terms in Aw The final second order contribution 
to the cross section through to O {!) is given by 

'^'^ - (^) {Aw + Asy + Ass) 



diTo \n 

'aV' 



TT/ 



+ L[Aioe'{^) + Bioe[^)+c 

The single logarithmic coefficients A, B and C are given by, 



(4.5) 



X 



A = 161og( ^— ^ ) - 16, (4.6) 



S = 8 Li2(l - x) - 8 Li2(x) + 12 log ( — + 4/(x) - 28 - |7^^ (4.7) 



,2 

l-xj ' ^■""^ 3'' 
93 5 

C - 6Li2(l - x) - 6Li2(x) + 2>f{x) + 6(3 - y - ^-k^ (4.8) 

Coefficients A and B agree with those obtained by Arbuzov et al. [9]^, while C is the 
main new result of the present work. 



^up to a slight misprint in the expression for B, denoted there by z\ 
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5. Summary 



In this form Eq. (4.5) cannot be used to compare directly with experiment. The 
reason is that an experimental set-up involves a complicated detector that is not 
represented by the simple energy cuts on the photons which we used here. Such 
effects have to be modeled by a Monte-Carlo calculation. However any uncertainty 
in the cross section is now reduced to O {a^) terms without enhanced logarithms. 
Given the final simplicity of the method one can wonder whether it is possible to 
use the results of [f 0] to extract the finite parts of the cross section. At present this 
appears to be not simply possible, as we have no way to tell whether the 1/e poles 
contained in H^'^^ (Eq. (2.5)) correspond to L or for instance L — 1. Presumably it is 
possible to establish this connection by further expanding the vertex functions. The 
alternative would be to calculate the box graphs within the mass regulator scheme. 
However this is quite a challenge. For application to analogous processes in QCD 
we remark however that this last uncertainty plays no role as it gets absorbed in the 
definition of the structure functions. 

The calculation presented in this paper removes a major obstacle to improvement 
of the precision of Monte Carlo programs for large angle Bhabha scattering at centre 
of mass energies of a few GeV [14] . The result can also be used at higher energies, 
provided that additional contributions due to Z-exchange diagrams are included. 
Whether it is possible to obtain them by an extension of the method we have used 
here is a question that requires further investigation. 
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